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Abstract This study was undertaken to estimate the effect 
of openings between cell walls on combined bound water 
and water vapor diffusion in wood. Using a newly devel¬ 
oped model, the radial and tangential moisture diffusion 
coefficients can be predicted depending on the opening 
area. The new model explicitly involves a term for water 
vapor diffusion through the openings, as well as a term for 
the combined diffusion of bound water and water vapor. A 
classical model developed by Stamm and Choong had 
higher longitudinal moisture diffusion coefficients than that 
in the parallel model at higher moisture content, which is 
inconsistent with the Wiener bound rule. The new model 
suggested in this article is useful for analyzing the experi¬ 
mental results and understanding the variability of the dif¬ 
fusion coefficients. 

Key words Bound water • Water vapor • Conductivity • 
Moisture diffusivity • Wiener bound 


Introduction 

Moisture movement below the fiber saturation point that 
involves bound water and water vapor is not an advanced 
topic in wood physics but a classical one. The integrated 
theoretical approach was initiated by Stamm 1 on the assump¬ 
tion that the moisture movement was divided into bound 
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water and water vapor and that each diffusion coefficient 
could be estimated. Therefore, the combined bound water 
and water vapor diffusion depends on the ratio of the cell 
wall to cavity volume and the geometrical structure of the 
cavity. This approach can simplify the understanding of the 
moisture movement in wood and it is still useful to wood 
scientists. Choong 2 extended the theoretical methods to the 
diffusion of moisture in the longitudinal direction. In this 
article, the model used by Stamm 1 and Choong 2 is referred 
to as the classical model. Siau 3,4 simplified the model by 
neglecting water vapor movement through the openings. 
The classical model and the Siau model cannot separate 
moisture diffusion coefficients in the radial and tangential 
directions. 

Water vapor flows through the pit or opening simultane¬ 
ously by convection and diffusion. The various models men¬ 
tioned above have some limitations in this respect. Stanish 
et al. 5 and Perre et al. 6 considered the convective and diffu¬ 
sive water vapor movement through openings for softwood 
drying models. There is great difference in the effective 
areas of the openings between the two models; that is, 
0.003-0.015 for Stanish et al. 5 and 0.001 for Perre et al. 6 Yeo 
et al. 7 calculated the internal and external resistances for 
moisture movement in wood and surface evaporation 
from it, and experimentally determined different diffusion 
coefficients in each wood direction: longitudinal, radial, and 
tangential. 

The accuracy of the simulation depends on the material 
parameters used. The apparent diffusivity by the experi¬ 
ments and the input value for simulation should have a solid 
physical meaning. Adapting the accurate diffusion coeffi¬ 
cients for the simulation, however, remains a serious 
problem. This is largely attributed to the natural variation 
within and between woods. In addition, it may be responsi¬ 
ble for the poor correlation between experiment and theory 
in some cases. Therefore, the objective of this study was to 
investigate the validity of the classical models and to develop 
a new model for combined bound water and water vapor 
diffusion. The transverse moisture diffusion coefficients for 
classical and Siau models are divided into radial and tan¬ 
gential components. 
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Modeling of combined diffusion coefficients 


The mass transfer equation is expressed by Eqs. 1 and 2. 


3(p m ra) 


dt 
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where p m and p™ are the wood density and cell wall density, 
respectively, based on oven-dry weight and moist volume 
(kg/m 3 ), m is the fractional moisture content, j h and / v are 
fluxes of bound water and water vapor (kg/m 2 s), and 0 is 
the porosity or void fraction. 

The porosity of wood may be calculated from wood and 
cell wall density using Eq. 3, assuming that the density of 
bound water is equal to the water density, the oven-dry cell 
wall density is 1530 kg/m 3 , and the pore volume remains 
constant with the moisture changes. 4 


A = a 2 = 1 P m = i Pq 1 + 1.53m 
^ pf 1 + mp 0 /1000 1530 


Vm<m fsp (3) 


where p 0 is the oven-dry wood density based on oven-dry 
weight and oven-dry volume (kg/m 3 ). 

It should be noted that porosity changes with moisture 
content. Bound water fluxes under isothermal conditions 
can be expressed by Fick’s law: 


7b = ~K h Vm = -pfDj^m 



where K b is the conductivity of bound water in cell wall of 
wood (kg/ms), and D B is the diffusion coefficient of bound 
water in the cell wall of wood (m 2 /s). 

For D b in Eq. 4, Stamm 8 first measured it in the longitu¬ 
dinal direction (T> BL ) after sealing off the lumens with 
molten bismuth. Siau 3 converted it to the diffusion coeffi¬ 
cient of the cell wall in the transverse direction (D BT ) by 
assuming D BT = D BL /2.5, and obtained the Arrhenius-type 
equation with the least-square method: 

£> bt = 7 x 10“ 6 exp [-(38500 - 29000m) / RT] (5) 

where R is the universal gas constant (8.314 J/molK), and T 
is temperature (K). The above diffusion coefficients would 
not represent all wood species and some researchers use 
other expressions. However, the equation has been used 
widely in the field of wood research. 9 

For water vapor fluxes, 


7v = ~K Y Vm = 


T) va M v p vs dh 


Vm 
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T RT dm 

where K w is the conductivity of vapor in the lumen of wood 
(kg/ms), D va is the bulk binary diffusivity of water vapor in 
air (m 2 /s), Tis the tortuosity (>1), h is relative humidity, M v 
is the molecular weight of water vapor (0.018kg/mol), and 
p vs is the saturated vapor pressure (Pa). 

In Eq. 6, the bulk binary diffusivity of water vapor in air 
by Dushman 4 is estimated by 

1.75 


D va = 2.2 x 10- 5 
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f T \ 
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I 273 J 
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where P is the total pressure of air and water vapor (Pa). 

There is little difference in cell wall density among wood 
species. The diffusion coefficient of water vapor in air and 
saturation vapor pressure can be determined using the rel¬ 
evant equations. Therefore, the combined movement of 
bound water with water vapor could be estimated if the 
information on the bound water diffusion coefficient, sorp¬ 
tion isotherm, and tortuosity are available. Tortuosity is 
defined as the ratio of the length of a true flow path for a 
fluid and the straight line distance and is always larger than 
1 . 0 . 

In this study, the tortuosity effect was neglected (t= 1). 
Thus, water vapor diffusion coefficient of air in the lumens 
based on the concentration of bound water in the cell wall, 
T) vc , is calculated from 


^va MyPy s dh 

p c m RT dm 



Neglecting the dimensional change with moisture 
content, the mass transfer equation, Eq. 1, can be simplified 
into 


d ^ = V(D h Vm) ( 9 ) 

dt 

The apparent directional moisture diffusion coefficients 
based on the concentration of bound water in wood ( D h ) 
can be calculated by dividing conductivities (K b ) by p m or 
p™( 1 - 0). Figure 1 shows primary paths for tangential and 
radial moisture transport through a cell wall (Fig. la) and 
secondary paths through a pit (Fig. lb). 

Electrical analog models for the paths have been devel¬ 
oped. Stamm 1 proposed an electrical model to simulate the 
moisture transfer in the hygroscopic range shown in Fig. 2a. 
Figure 2b is a simplified model that integrates the secondary 
paths through a pit. Because the effective area of the pit 
openings in softwoods is small and the resistance is large, 
Siau 4 neglected their effects and adopted the model shown 
in Fig. 2c. The contribution of the pits in hardwood is 
expected to be less than that in softwood because of smaller 
pit openings. As Siau stated, however, the effect of pit open¬ 
ings may be important in woods with high permeability and 
high specific gravity, particularly at low moisture contents 
and high temperature. 
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Fig. la, b. Primary paths for tangential and radial moisture transport 
through a a cell wall, and b secondary paths through a pit 
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For the classical model as shown in Fig. 2b, the resistance 
and conductance of the parallel-series model are 



gb 
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For the Siau model of Fig. 2c, the flux direction is assumed 
to be perpendicular to that of the classical model/ It should 
be noted that the area of cross walls and the length of side 
walls is different from the classical model. In this study, the 
geometry of cell walls used in the classical model is adopted. 
However, the path length through the pit for the new model 
is 1.0 but 1 - a for the classical model. Therefore, the electri¬ 
cal analog for the paths is somewhat different from the 
models. 

Similarly, for the new model of Fig. 2d, 
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Therefore, conductivity is defined by 


( 12 ) 

(13) 


K b gb 




(a) 


92 

AAAr 


94 

A/W 


9i 

A/W 


96 

95 rAAAn 

AVA 97 a j 
AW 


(b) 


94 

AAAr 


92 

AW 


9i 

AMn 

93 

WAV 


(c) 


92 

9i rAAAn 


AAAr 


94 


(d) 


94 

AAAr 

93 

AAAr 


92 gi 

AM—VVW 


Fig. 2a-d. The electrical analog models for the paths, a Stamm model, 
b Stamm’s simplified model, c Siau model, d new model. Conductance: 
g 1? cross wall; g 2 , lumen; g 3 , pit or ray; g 4 , side wall; g 5 , pit chamber; g 6 , 
pit membrane wall; g 7 , pit pore 


where A e is cross-sectional area perpendicular to flux direc¬ 
tion and L e is length in the flux direction. 

For a mixture of solid A and void B as shown in Fig. 3, 
the diffusion coefficients can be derived using the mixing 
rules. The series model for conductance is obtained by the 
harmonic mean. 


1 _1-0 , 0 
gs gA g B 



The parallel model is calculated using the weighted 
mean. 


gp=(l-0)gA + 0gB (16) 

For a dispersed mixture, it may be calculated by the 
Maxwell model. 10,11 


gm gB 


" 2g B +gA-2(i-<ft)(g B -gAr 

2 g B + g A + (l- 0 )(gB-gA) J 



The series and parallel model provide the so-called Wiener 
bounds, which define the upper and lower bounds for the 
effective diffusion coefficient of a mixture. The real diffu¬ 
sion coefficient of any mixture should be between these 
boundaries. Figure 3d-f shows schematic representations of 
the three directional diffusion models dealt with in this 
study. 


Results and discussion 

The conductivities of wood in the principal directions can 
be calculated by substituting the conductivities of the cell 
components (Table 1) into Eqs. 11, 13, and 14. The results 
are as shown in Table 2. 

Tangential and radial models are proposed, with a unit 
cube of material measured between two parallel surfaces 
(A e = 1.0 and L e = 1.0). A longitudinal model with a di¬ 
mension of unit area and length of L + 1 - a is proposed 
(A e = 1.0 and L e = L + 1 - a). For the classical model (the 
Stamm and Choong models), water vapor diffusion into the 
cavity and through openings is always involved in a com¬ 
bined term with bound water. Therefore, the classical model 
cannot explain the effect of water vapor through the open¬ 
ings in a multiphase moisture transport model. 5,6 However, 
the new model explicitly involves a term for water vapor 
diffusion through the openings as well as the diffusion of 
combined bound water and water vapor. 


Table 1. Conductance of cell components in the principal directions 


Direction 

Cross walls (g,) 

Lumen (g 2 ) 

Pit or ray (g 3 ) 

Side walls (g 4 ) 

Tangential 

Pc iAt A — p ) 

1 -a 

K v (a-p) 

a 

mat 

vl -a) 

Pc DbtO- ~ a ) 

Radial 

Pc D bt u 

1 — a 

K 

K v r 

Pc DbtO- ~ a ~ r ) 

Longitudinal 

pT D BL a (a - p) 

K v a(a-p) 

K w pa ( K v pa ) a 

Pc 1 A) 


1 -a 

L 

L+l-a\ 1 —aJ 

L + l — a 


a Classical model 
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Fig. 3a-f. Schematic 
representation of the diffusion 
models for a mixture and wood, 
a Series model, b parallel model, 
c dispersed mixture model, 
d tangential model, e radial 
model, f longitudinal model, p , 
Ratio of an effective pit opening 
area to a cross-wall or end-wall 
area; a , ratio of an effective 
diameter of lumen to diameter 
of cell; r, ratio of an effective 
ray area to a cell wall area; L, 
ratio of fiber length to diameter 
of cell 
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Table 2. Conductivity of wood in the principal directions 
Model Conductivities 


New 

Tangential K 

Radial K 

Longitudinal K 

Classical 

Tangential K 

Longitudinal K 


bT — 


bR - 


bL — 


bT — 


bL — 


C a~p)p?D 


BT 


1 - a + ap™D BT / D v 
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(a - p) [(a - p) p?P BT + p Ac ] 
a-p-a 2 + lap + a(a-p)p™D m /D 

(a - p) (L +1 + a) [(a - p) p c m D BL + pD, c ] 

1 -a + p(l-lla + L/a) + L(a- p)p™D BT / D y 


BT 


+ ( 1 - 


a 



In extreme cases where p-a , the diffusion model must 
be a parallel model because the area of the cross wall is 
close to zero. For the new model, 

K bT = p?D BT (l-a) + D vc a (18) 

^bL = D vc a 2 + p™D BL (l- a 2 ^) (19) 

In the classical model for p = a, 

= Pc 1 £*Bt(1 _ a ) (^0) 


K bL = p?D BL (l-a 2 ) (21) 

Figure 4 shows that tangential and longitudinal moisture 
diffusion coefficients with oven-dry specific gravity G 0 = 0.8 
and p = a at 30°C, which result from Eqs. 18-21. The 
new model’s tangential diffusion coefficient (£> bx ) is nearly 
equal to that of the parallel model, and the longitudinal 
diffusion coefficient (D bL ) is lower than that of the parallel 
model due to the difference in the estimated porosity. 
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Moisture content 

Fig. 4. Tangential and longitudinal moisture diffusion coefficients with 
G 0 = 0.8 and p = a at 30°C. Solid line, new model; dotted line, classical 
model 



Moisture content 

Fig. 5. Tangential and longitudinal moisture diffusion coefficients with 
G 0 = 0.8 and p = 0 at 30°C. Solid line, new model; dotted line, classical 
model; dash dotted line, Siau model 


However, D hL and D hT of the classical model are closer 
to those of the series model, which is erroneous and 
unrealistic. 

The conductivities by the new model for p = 0 are 




ap™D 


BT 


1 — a + p c D BB a / Z) vc 

a p m D BT 


+ (l-a)p c D BT 
p m D BT 


+ 


1 — a 1 — a-\- ap c D B ^ / £>vc 1T Cl 


a 


(L + l-a)p m D BL 


1 -a A 1 - a + Lp^/^BL ID 


+ p m D 


BL 


VC 




For the classical model for p - 0, the conductivities are 


a _ p m D bt _ | P m D BT 

1-a 2 l-a + ap^G^BT/^vc 1 + 0 

a 2 (L + l-a)p m Z) BL ^ mr> 

~ ; sr; 77^~ + p ^bl 

1 - a 1 - a + Lap c T) B L /Dye 


(24) 

(25) 


The tangential conductivity of the classical model is the 
same as that of the new model, but the longitudinal conduc¬ 
tivity is slightly different between the classical and new 
models. 

Figure 5 shows D hT and D hL for p - 0. There is little dif¬ 
ference between the Siau model and the new model in 
terms of D hh . The D hL in the classical model is always higher 
than that in the new model and the Siau model. The classi¬ 
cal model has a higher D hL than that in the parallel model 
at higher moisture content, which is inconsistent with the 
Wiener bound rule. The D hT or transverse moisture diffu¬ 
sion coefficients in the Siau model are always 1/a times 
higher than those in the new model because of the larger 
area of the cross walls and the smaller area of the side walls. 
As the wood density increases, therefore, the difference 
between the Siau model and the new model increases. 

Figures 6 and 7 show the moisture diffusion coefficients 
at various p, respectively, for G 0 = 0.4 G 0 = 0.8. The effect 
of the openings is more important in woods with higher 



Moisture content 


Fig. 6. Moisture diffusion coefficients with G 0 = 0.4 at various p. Solid 
line,p (or r) = 0; dotted line,p (or r ) = 0.001; dash dotted line,p (or r) 
= 0.01 



Moisture content 

Fig. 7. Moisture diffusion coefficients with G 0 = 0.8 at various p. Solid 
line,p (or r) = 0; dotted line; p (or r ) = 0.001; dash dotted line,p (or r) 
= 0.01 










348 


permeability and higher specific gravity, particularly at 
lower moisture contents and higher temperature. For woods 
with higher permeability, D bT and D bR are nearly constant 
with moisture content. Depending on the magnitude of p 
and r, the ratio of radial and tangential diffusion coefficients 
may change. Because the tangential model and the radial 
model have the same schematics when p and r - 0, the 
marks of D bT and D bR in the figures overlap each other. 


Conclusions 

The flux and path length of water vapor were corrected in 
terms of the combined bound water and water vapor diffu¬ 
sion coefficients that were originally deduced by Stamm. 1 
The transverse moisture diffusion coefficients in the Siau 
model were always higher than those in the new model by 
a factor equivalent to the inverse of the ratio of effective 
diameter of the lumen to the diameter of the cell, because 
of the larger area of the cross walls and the smaller area 
of the side walls. The new model can predict radial and 
tangential diffusion coefficients. The new model explicitly 
involves a term for water vapor diffusion through the open¬ 
ings as well as combined diffusion of bound water with 
water vapor. It may be helpful for analyzing the experimen¬ 
tal results of the diffusion coefficients and understanding 
the variability. In addition, it would be useful for calculating 
the effect of moisture flux on heat transfer. For woods 
with larger openings or higher permeability, the apparent 
diffusion coefficient might be near constant with changing 
moisture content. 
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